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Abstract. In this work, we studied the stehar and baryonic Tuhy-Fisher 
relations by using hydrodynamical simulations in a cosmological frame- 
work. We found that supernova feedback plays an important role on 
shaping the stellar Tully-Fisher relation causing a steepening of its slope 
at the low-mass end, consistently with observations. The bend of the re- 
lation occurs at a characteristic velocity of approximately 100 km s~^, in 
concordance with previous observational and theoretical findings. With 
respect to the baryonic Tully-Fisher relation, the model predicts a linear 
trend at z ~ with a weaker tendency for a bend at higher redshifts. 
In our simulations, this behaviour is a consequence of the more efficient 
action of supernova feedback at regulating the star formation process in 
smaller galaxies. 

Resumen. En este trabajo, estudiamos las relaciones de Tully-Fisher 
estelar y barionica mediante la utilizacion de simulaciones hidrodinamicas 
dentro de un marco cosmologico. Encontramos que la retroalimentacion 
al medio por supernovas juega un rol importante en la determinacion de 
la forma de la relacion de Tully-Fisher estelar incrementando su pendien- 
te en el sector menos masivo, consistentemente con las observaciones. El 
cambio de pendiente de la relacion ocurre para una velocidad caracte- 
ristica de aproximadamente 100 km s~^, en concordancia con hallazgos 
observacionales y teoricos previos. Con respecto a la relacion de Tully- 
Fisher barionica, el modelo predice un comportamineto lineal a z ~ 
con una tendencia debil a curvarse hacia corrimientos al rojo altos. En 
nuestras simulaciones, este comportamiento es una consecuencia de la ac- 
cion mas eficiente de la retroalimentacion al medio por supernovas en la 
regulacion de la formacion estelar de las galaxias mas pequehas. 



1. Introduction 

The correlation between the luminosity and the rotation velocity for spiral galax- 
ies (TFR, Tully &; Fisher 1977) has been widely studied during the last decades. 
However, it is now accepted that the TFR is a result of the more fundamental 
scaling relation between the rotational velocity of a galaxy and its stellar mass 
(sTFR) or baryonic mass (bTFR). 
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In particular, there is evidence that the sTFR exhibits a change in its slope 
at around ~ 90 km s^^ in the sense that slow-rotators have lower stellar masses 
than those derived from the extrapolation of the linear fit for fast-rotators (Mc- 
Gaugh et al. 2000; Amorm et al. 2009). Moreover, McGaugh et al. 2010 have 
recently reported that there is also a bend in the bTFR but located at a lower 
rotation velocity. 

In this work, we performed hydrodynamical simulation in a ACDM universe 
to study the role of supernova (SN) feedback on the origin of the shape of the 
sTFR and bTFR and on their evolution with redshift. 

2. Simulations and sample selection 

We performed numerical simulations consistent with the concordance ACDM 
universe with = 0.3, Oa = 0.7, Qb = 0.04 and Hq = 100 /i"^ kms'^ Mpc"^ 
with h = 0.7. A version of the chemical code GADGET-3 including treatments 
for metal-dependent radiative cooling, stochastic star formation, a multiphase 
model for the ISM and SN feedback (Scannapieco et al. 2006) was employed. 
The simulated volume corresponds to a cubic box of a comoving 10 Mpc 
side length. The simulation has a mass resolution of 5.93 x IO^Mq/i"^ and 
9.12 X 10^Mq/i~^ for the dark and gas phase, respectively. 

Simulated disc-like galaxies were identified following the methods describe 
by De Rossi et al. (2010). The mean properties of galactic systems were esti- 
mated at the baryonic radius (-Rbar)) defined as the one which encloses 83 per 
cent of the baryonic mass of the systems. We found that the tangential velocity 
of these systems constitutes a good representation of their potential well so that, 
for the sake of simplicity, we used the circular velocity estimated at i?bar as the 
kinematical indicator for our study. 

3. Results and discussion 

With respect to the local sTFR, fast-rotators (100 km s~^ < V < 250 km s~^) de- 
scribe a linear relation of the form log(M*/M0/i-^) = (3.68±0.09) log{V/ 100 km s"^) 
+(9.42 ±0.26), which is in general good agreement with observations (e.g. Mc- 
Gaugh et al. 2000). Moreover, our model is also able to reproduce the observed 
steepening of the relation in the case of slow-rotators {V < 100 km s~^; see 
De Rossi et al. 2010 for details), in agreement with previous observational and 
theoretical works (Larson 1974; Dekel & Silk 1986; McGaugh et al. 2000). In 
order to analyse at what extent SN feedback might be responsible for the bend 
of the sTFR, we suppress this mechanism from our simulations. We found that, 
when SN feedback is turned-off, a single linear relation is recovered with a flatter 
slope (~ 2.97), suggesting the important role of SNs in the determination of the 
shape of the sTFR. 

In the case of the local bTFR, we obtained a linear trend for the whole 
range of circular velocities resolved by these simulations (40 km s~^ < 1^ < 
250 km s~^). This result does not disagree with the findings of McGaugh et 
al. (2010) since they reported that the bend of the bTFR is at a lower rotation 
velocity {V ~ 20 km s~^). Furthermore, the simulated bTFR follows a relation 
of the form log{Mb/MQh-^) = (3.23 ± 0.08) log{V/ 100 km s^^) + (9.56 ± 0.23), 
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consistently with observations for late and early type galaxies (De Rijcke et al. 
2007; Gurovich et al. 2010). 

By analysing the simulated sTFR and bTFR within the redshift range < 
z < 3, we found that neither of them show a significant evolution of the slope 
at the high-mass end. However, in the same redshift range, we detected an 
increase of the zeropoint by ~ 0.44 dex and ~ 0.33 dex for the sTFR and bTFR, 
respectively. 

In order to explore how these trends arise in our SN model which does not 
introduce scale-dependent parameters, we analysed its effects on the gas-phase of 
simulated galaxies. Firstly, we investigated the importance of SN-driven outflows 
by calculating the fraction fi = ilmAfj/il{,Mvir for each galaxy at z = 0, where 
Mvir is the total mass within the virial radius and Mj denotes the stellar (/*) 
or the baryonic (/;,) component within -Rbar- Our results indicate that is an 
increasing function of the circular velocity suggesting that SN feedback might 
be strongly regulating the star formation efficiency (eSFR) of smaller galaxies. 
To analyse this issue, we compared the correlation between the eSFR and the 
circular velocity of the systems in the case of the SN-feedback model with the 
one obtained from a feedback-free run. As expected, the SN feedback model 
predicts lower eSFRs at a given velocity with the stronger differences at the 
low-mass end of the relation, consistently with the steepening of the sTFR for 
slow-rotators. By studying the fraction fh as a function of the circular velocity, 
we found a weaker correlation than in the case of which accounts for the 
trend to recover a single slope for the bTFR. However, the fact that fh < 0.6 
for the whole sample suggests that galactic winds are important over the entire 
range of simulated stellar masses. Finally, we estimated the fraction f^^^ defined 
as the ratio between the simulated baryonic mass within the virial radius and 
the expected one. By comparing and f™, we determined that an important 
amount of the missing baryons within galaxies can be found in the surrounding 
halo (a detailed discussion can be found in De Rossi et al. 2010). However, 
given that /^"' is within the range 0.1 — 1 for the whole sample, it is clear 
that a significant fraction of the gas-phase of simulated galaxies is blown away 
as a consequence of very efficient galactic winds. And, these effects are more 
prominent in smaller systems. 

Basically, in our model, star formation is triggered mainly when the gas 
gets cold and dense ^. As a consequence of star formation, SN energy is released 
heating up the surrounding cold ISM. The smaller virial temperatures of slow 
rotators lead to a more efficient transition of the gas from the cold to the hot 
phase, and hence, generates a decrease of the star formation activity in these 
systems. However, the cooling times for these galaxies are still too short com- 
pared to the dynamical times to allow this hot phase to be stable and the gas 
can return to the cold phase in short time-scales. In fact, SN feedback leads to 
a self-regulated cycle of heating and cooling generating an important regulation 
of the star formation process in smaller galaxies. In the case of fast rotators, the 



^In this simulations, the cold phase is defined as the gas component with temperature T < Tc 
where Tc — 8 x 10'' K and density p > O.lpc where pc is 7 x 10~^^g cm~^. Otherwise, the gas 
is classified as hot phase. The reader is referred to Scannapieco et al. 2006 for more details 
about the model. 
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hot phase is estabhshed at higher temperature due to the continuous injection 
of SN energy and the increase of the virial temperature of the dark halos hosting 
them. And, since in these systems the coohng times get longer compared to the 
dynamical times, the hot gas is able to remain in this phase. However, at this 
stage, galactic winds will be more difficult to be triggered since the SN energy 
accumulated by the cold phase is not enough to match that of its nearby hot 
environment, which is a condition required by the model to promote particles 
to the hot phase. Meanwhile, this gas remains available for star formation and, 
consequently, SN feedback is not efficient at regulating the star formation in 
larger galaxies. Interestingly, this model predicts that the transition from effi- 
cient to inefficient gas cooling occurs at around the same characteristic velocity 
of 100 km where the sTFR bends and is also in agreement with previous 
observational and theoretical findings. 

4. Conclusions 

We have studied the sTFR and bTFR by using cosmological simulations in a 
ACDM universe. Our results suggest that SN feedback seems to be crucial to 
reproduced the observed bend of the sTFR as a consequence of the more efficient 
action of SN feedback in the regulation of the star formation activity in smaller 
galaxies. Without introducing scale-dependent parameters, the model predicts 
that the bend occurs at a characteristic velocity of ~ 100 km s~^, consistently 
with previous observational and theoretical works. The reader is referred to De 
Rossi et al. (2010) for more details about this work. 
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